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Abstract 

We explore the cosmological evolution in the exponential gravity f{R) = i? + ci (l — e~'^^^) 
(^ i (ci 2 = constant). We summarize various viability conditions and explicitly demonstrate that the 

f-^ '_ late-time cosmic acceleration following the matter-dominated stage can be realized. We also study 

(D ■ 
^/) • the equation of state for dark energy and confirm that the crossing of the phantom divide from the 



a^ 



phantom phase to the non-phantom (quintessence) one can occur. Furthermore, we illustrate that 



r~N I the cosmological horizon entropy globally increases with time. 
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I. INTRODUCTION 



There exist two representative approaches to account for the current accelerated expan- 
sion [l|-l3| of the universe |#^9|. One is to introduce "dark energy" in the framework of 
general relativity. The other is to consider a modified gravitational theory, such as f{R) 
gravity. 



Several viable theories of f{R) gravity have been constructed; e.g., power-law |10| . Ill| . 
Nojiri-Odintsov [l2|, Hu-Sawicki 13[, Starobinsky's [lj|, Appleby-Battye 15[, and Tsu- 



jikawa's [16| models (for more detailed references, see a recent review on f{R) gravity [7[). 



0; 



It is known that these models can satisfy the following conditions for the viability: (i) positiv- 
ity of the effective gravitational coupling, (ii) stability of cosmological perturbations |l7H20l|. 



(iii) asymptotic behavior to the standard A-Cold-Dark-Matter (ACDM^ 



10 



21 



model in the large 



22| . (v) constraints 



curvature regime, (iv) stability of the late-time de Sitter point 
from the equivalence principle, and (vi) solar-system constraints 23 1. 

Recently, an interesting model of f(R) = -R + Ci (l — e~^^^), called "exponential grav- 
ity", has been proposed in Refs. 24|-|26| with ci^2 being constants. The important feature 
of the exponential gravity is that it has only one more parameter than the ACDM model. 
The constraints from the violation of the equivalence principle |27| and cosmological ob- 



servations |28[ on the exponential gravity have been examined. The exponential gravity 
in the framework of f(R) gravity has been extended to a gravitational theory in terms of 
the torsion scalar 29[ (for a related work on torsion gravity, see 30[). We note that the 
cosmological dynamics in the gravitational theory consisting only of the exponential term 
without the Einstein- Hilbert one has also been studied in Ref. 31|. 



In this paper, we explicitly investigate the cosmologica 
gravity model given by Cognola et al. 25| and hinder 



evolution in the exponential 
in more detail by using the 



analysis method in Ref. [13| . We also check the above six viability conditions for the model. 
In particular, we demonstrate that after the matter-dominated stage, the current accelerated 
expansion of the universe and the crossing of the phantom divide from the phantom phase 
to the non-phantom (quintessence) one can be realized. It is interesting to note that the 
crossing of the phantom divide is implied by the cosmological observational data 32| , while 
the exponential gravity is a ghost free theory. In addition, we illustrate that the cosmological 
horizon entropy globally increases with time. We use units oi k-Q = c = h = 1 and denote 



the gravitational constant SnG by k."^ = Sn/Mpi^ with the Planck mass of Mpi = G ^/^ = 
1.2 X 10^9 GeV. 

The paper is org anized as follows. In Sec. II, we review the model of the exponential 



gravity in Refs. 25|, |26| and summarize its viability conditions. In Sec. Ill, we explore the 
cosmological evolution of the model. We examine the horizon entropy in Sec. IV. Finally, 
conclusions are given in Sec. V. 

II. EXPONENTIAL GRAVITY 

A. The model 



The action of f{R) gravity with matter is given by 

:/CR) 
2k' 



d X\/ —g ^2 1" -< matter (fi'^j/ 5 Tmattcr) ; (2-1) 



where g is the determinant of the metric tensor g^^^ /matter is the action of matter which is 
assumed to be minimally coupled to gravity, i.e., the action / is written in the Jordan frame, 
and Tmatter deuotcs matter fields. Here, we use the standard metric formalism. 

Taking the variation of the action in Eq. (12.11) with respect to g^^y^ one obtains [6| 

FG,, = /.^T^-"-) - 1^^, {FR - /) + V^V.F - g.^DF , (2.2) 

where G^^ = R^u — (1/2) QfiuR is the Einstein tensor, F{R) = df{R)/dR, V^ is the covariant 
derivative operator associated with g^^^, D = g^'^V ^Vu is the covariant d'Alembertian for a 
scalar field, and T^™^ '^^ is the contribution to the energy-momentum tensor from all perfect 
fluids of generic matter. 

In this paper, we concentrate on the exponential gravity in Refs. 25|, |26| , given by 

/(/?) = i?-/3i?s(l-e-^/^=) , (2.3) 

where ci = — /3-Rs and C2 = R^^- Note that R^ corresponds to the characteristic curvature 
modification scale. 

B. Viability conditions on exponential gravity 

For the model of the exponential gravity in Eq. (12. 3p . it is straightforward to show that 
the conditions for the viability can be satisfied, which are summarized as follows: 



(i) Positivity of the effective gravitational coupling 

When /3 < e^/^% F{R) = 1 — I3e~^/^^ > 0. This is required for the positivity of the 
effective gravitational coupling Gcs = G/F{R) > to avoid anti-gravity. In the sense of the 
quantum theory, the graviton is not a ghost. 
(ii) Stability of cosmological perturbations 

When /3 > and i?s > 0, f"{R) = F'{R) = {(3/R,) e-^/^ > 0, where the prime 
denotes differentiation with respect to R. This is required for the stability of cosmological 



perturbations 18|-|20| . In the sense of the quantum theory, the scalaron, which is a new 



scalar degree of freedom in f{R) gravity, is not a tachyon 14 1. 

(iii) Asymptotic behavior to the ACDM model in the large curvature regime 

Since f{R) — /?—)■ —f3Rs = constant for R/Rs ^ 1, this model is reduced to the ACDM 
model in the large curvature regime R/Rs ^ 1. Such a behavior is necessary for the presence 
of the matter-dominated stage. 
(iv) Stability of the late-time de Sitter point 

When /3 > 1, < m{R = Rd) < I m, where m = Rf"{R)/f'{R) = RF'{R)/F{R) 
and -Rd = 2/(-Rd)/-^(-Rd) is the value of the scalar curvature at the de Sitter point. This 



condition is required for the stability of the late-time de Sitter point |10l . |21| . |22| . The 
quantity m characterizes the deviation from the ACDM model because ttt, = for the 
ACDM model. In the exponential gravity, by using m{R) = /3 (R/Rg) e~^/^^/ (l — /3e"^/^=) 
and (3 = (Rd/Rs) / [2 - (2 + Rd/Rs) e~^-^/^] , one finds that m{Rd) < 1 for Rd/Rs > 0. 
Hence, m{Rd) < 1 for /3 > 1. 
(v) Constraints from the violation of the equivalence principle 

It is known that f{R) gravity can satisfy local gravity constraints from the violation of 



the equivalence principle under the chameleon mechanism 33|, |3J]. By making the following 



conformal transformation [35|: g^^ — )■ g^^ = 'E?g^^, the action of f{R) gravity in Eq. (12. ip 
can be rewritten in the Einstein frame, where S^ = F = eV^/^'"?^ with the scalar field 
0. In what follows, a tilde represents the quantity in the Einstein frame. We consider a 
spherically symmetric body with radius fc in the Minkowski space-time. Here, f is the 
distance of the center of the body, and p* = e~v '^'^"'^p is a conserved matter density in 
the Einstein frame with p the energy density of matter in the Jordan frame. We assume 
that a spherically symmetric body has constant densities of p* = p^ and Pout(^ Pin) inside 



(f < fc) and outside (f > fc), respectively. In this case, the effective potential has two 
minima at the field values 0in and 0out satisfying the conditions dVea{(f)in)/d(t> = and 
dVcs {(f> out) / d(t) = with a heavier mass squared mf^ = d'^Vcs{(l)in)/d(f>'^ and a lighter mass 
squared m^^^ = d'^Vcf[{(pout}Jd(j)'^ , respectively. The thin-shell parameter is defined as eth = 
— /^ (0out ^ 0m) / ("\/6$c) 3^, where $c = GMc/fc is the gravitational potential at the 



surface of the body and Mc = (47r/3) r^Pm- 

The tightest experimental bound on eth obtained from the violation of the equivalence 
principle for the accelerations of the Earth and the moon toward the Sun is given by eth, ® < 



2.2x10-6 



33, 



38 



This is the thin-shell parameter for the Earth. By using the value of the 
gravitational potential for the Earth $c,® = 7.0 x lO^"*^" and |0out,®| ^ l^m,®!, the condition 

3l|. The field value 0out,® can be found by 



on eth,® is reduced to |fi;0out,©| < 3.7 x 10 ^^ 

solving dVcs {(pout) /d(t> = with p* = pout, which gives R ~ /tVout^ 

For the exponential gravity, Kcpout - -A/^/Je"'''''""*/^'' J27| and (3RjRo ^ 11^°^ 



m , 



where Rq ~ ^'^H^ is the current scalar curvature, Hq is the current Hubble parame- 
ter, fiin = Pm /Pcrit ^^ ^^^ currcut density parameter of non-relativistic matter (cold 
dark matter and baryon), pL is the energy density of non-relativistic matter at the 
present time, and pl.J^^ = 3Hq/k'^ is the critical density. As a consequence, by using 
Pcrit — IQ-^^g/cm^ and the homo gen eous baryon/dark matter density pout — 10-^^g/cm^, 



we find K</)out ~ —/3exp (—10^/3) 27|. When /3 > 1, which is the stability condition for 
the late-time de Sitter point in the exponential gravity, the above constraint on |fi;0out, ®| 
is satisfied very well. For example, if /3 = 1.1, |fi;<;/)out| = O{10~^^^^^). In what follows, the 
superscript (0) denotes the present value, 
(vi) Solar-system constraints 

The bound on the thin-shell parameter coming from the solar-system constraint eth,® < 
2.3 X 10"^ |7| is weaker than that from the violation of the equivalence principle eth,® < 
2.2 X 10"^ shown above. 

III. COSMOLOGICAL EVOLUTION 

We assume the fiat Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time with the 
metric, 

ds^ = -dt^ + a^{t)dx\ (3.1) 



where a{t) is the scale factor. From Eq. (12. 2p . we obtain the following gravitational field 
equations: 

3FH^ = k'^pm + ^{FR-f)- 3HF , (3.2) 

-2FH = K^ (pM + Pm) + F-HF, (3.3) 

where H = a/a is the Hubble parameter, the dot denotes the time derivative of d/dt, 
and pm and Pm are the energy density and pressure of all perfect fluids of generic matter, 
respectively. 

Equation (13.21) can be rewritten to 



H^-{F-l)\H^ + H'\ + ^-U-R) + H'F'^='^, (3.4) 

\ dma /D a ma 3 

while the scalar curvature R is expressed as 

R = 12H^ + Qh4^ . (3.5) 

a In a 



To solve Eqs. (13. 4p and (13. 5p . we introduce the following variables |13| : 
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with 
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(3.6) 
(3.7) 

(3.8) 
(3.9) 



.(0) 



where pde is the energy density of dark energy and pr is the energy density of radiation 
at the present time. In our analysis, the contribution from radiation is also taken into 
consideration. Equations (13. 4 p and (13. 5p are reduced to a coupled set of ordinary differential 
equations 

dyu Vr . /Q 1 n^ 

^-%H, (3.10) 



din a 3 

dVR ^„_3 1 1 



9a 



din a i/f/ + a '^ + xa ^m?F' 



X 



\ / 1 -'- -S —4 1 -'- / — ^ 

yH-[F -l)\ -Vr - yif - -a - xa 



6 2 /6m 



m^ 



(3.11) 



The equation of state for dark energy wbe = -Pde/pde is given by 

wbe = -1-- — -7] — , (3.12) 

derived by the continuity equation 

Pde + 3if (1 + m;de) Pde = . (3.13) 

On the other hand, the effective equation of state Wcs is defined as 

-..= -1-3^ = —, (3.14) 

where ptot = Pde + Pm + Pr and Ptot = -Pde + -Pm + -Pr are the total energy density and 
pressure of the universe, respectively. Here, Pde, -Pm (= 0) and Pj. are the pressure of dark 
energy, non-relativistic matter and radiation, respectively. 
Combining Eqs. f l3.10p and (13. lip , we obtain 

Ji-r, ^ J2yH + Js = , (3.15) 



d {hi. ay din a 



where 



1 2 — F 

yu + ci + XO om''r ' 

yn + a + X*^ 6m^F' 

In Figs. 1, 2 and 3, we depict the cosmological evolutions of the density parameters of dark 
energy 1],, . p^^/p^ non-relativistic matter Q^ ^ pjf/^ and radiation n^ ^ pjf/^ 
as functions of the redshift z = 1/a — 1 for /3 = 1.1, /3 = 1.8 and /3 = 2.5, respectively. 
In the high z regime [z > 3.0), the universe is at the matter-dominated stage (^2^ > ^de, 
Q^ ^ Qj.). As z decreases, the dark energy becomes dominant over matter for z < z^-e, 
where zqe is the crossover point in which Qj^e = ^m- Explicitly, we have z-[)e = 0.55, 0.47 
and 0.45 for /3 = 1.1, 1.8 and 2.5, respectively. The values of zde become smaller for the 
larger values of (3. At the present time {z = 0), (f]g^, f]S\ f^r°^) = (0.77,0.23,7.0 x 10"^), 
(0.76,0.24,7.3 x 10"^) and (0.75,0.25,7.3 x 10"^) for (3 = 1.1, 1.8 and 2.5, respectively. 
In Fig. 4, we also show the cosmological evolution of Q,- for /3 = 1.8. The qualitative 
behaviors of Q^ for P = 1.1 and 2.5 are similar to that for /3 = 1.8. Thus, the current 




FIG. 1: Cosmological evolutions of Ode (solid line), Om (dashed line) and Qj- (dashed and single- 
dotted line) as functions of the redshift z for (3 = 1.1. 



accelerated expansion of the universe following the matter-dominated stage can be realized 
in the exponential gravity. 

We note that in solving Eq. (13.181) numerically, we have taken the initial conditions at 
z = Zi as Hh {z = Zi) = 3.0 and dyn/dlna^z = z^) = 0, where Zi = 4.0, 3.5 and 3.0 for 
P = 1.1, f3 = 1.8 and (3 = 2.5, respectively. The values of Zi have been chosen so that 
RF'{z = Zi) ~ 10"^'^, i.e., the exponential gravity a.t z = Zi can be very close to the ACDM 
model, in which RF' = 0. Since R/Rs ^ 1 in the high z regime [z ~ zi), PRs/rh"^ ~ Quh- 
Consequently, the value of the combination (3Rs is set as (3Rs — 18/70^2^. Therefore, we 
have only one free parameter /3 in the exponential gravity in Eq. (12. 3p . Furthermore, from 
Eq. (I3.10p we see that i/r = 12yH at z = Zi and it follows from Eq. (I3.12p that t^DE = — 1 at 



8 




FIG. 2: Legend is the same as Fig. 1 but for (3 = 1.8. 



>(o) 



z = Zi. All numerical calculations have been executed for Qm' = 0.26 J2|. 

The cosmological evolution of the equation of state for dark energy wde in Eq. fl3.12p 
is shown in Fig. 5. From the figure, we see that wde starts at the phase of a cosmological 
constant Wde = — 1 and evolves from the phantom phase (wde < —1) to the non-phantom 
(quintessence) phase (wde > — !)• The crossing of the phantom divide occurs at z = Across, 
where -Across = 0.78, 0.57 and 0.46 for (3 = 1.1, 1.8 and 2.5, respectively. The values of 
Across become smaller for the larger values of /3. Moreover, the present values of wde are 
wbe{z = 0) = —0.85, —0.93 and —0.97 for /3 = 1.1, 1.8 and 2.5, respectively. Since (3Rs 
is a constant, the larger /3 is, the closer the exponential gravity is to the ACDM model. 
The results on wde are qualitatively the same as the analysis in Refs. 26|, |28|] . Thus, the 
crossing of the phantom divide from the phantom phase to the non-phantom one can be 
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FIG. 3: Legend is the same as Fig. 1 but for /3 = 2.5. 



realized in the exponential gravity. W e remark that the sijnilar behaviors can occur in 
Hu-Sawicki 



13 



42|, Appleby-Battye J43|, and Starobinsky's J44| models as well. 



In Fig. 6, we also illustrate the cosmological evolution of the effective equation of state 
Woff in Eq- fl3.14p . The present values of Wcg are Wf,<s{z = 0) = —0.65, —0.71 and —0.74 for 
/3 = 1.1, 1.8 and 2.5, respectively. We remark that w^s does not cross the line of the phantom 
divide unlike w;de due to the null energy condition ptot+-Ptot = PDE+Pm+Pr+-PDE+-Pm+-Pr > 
0. 

Finally, we mention that an /(-R) gravity model with realizing a crossing of the phantom 
divide from the non-phantom phase to the phantom one, which is the opposite transition 
from the above one, has been reconstructed in Ref. 45|]. In addition, the behavior of /(-R) 
gravity with realizing multiple crossings of the phantom divide 46| and that of /(-R) gravity 
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FIG. 4: Cosmological evolution of 0,^ (solid line) as a function of the redshift z for /3 = 1.8. 



around a crossing of the phantom divide by taking into account the presence of cold dark 
matter 47|] have also been explored. 



IV. HORIZON ENTROPY 



In Ref. 



48|, it is shown that it is possible to obtain a picture of equilibrium thermody- 



namics on the apparent horizon in the FLRW background for f{R) gravity as well as that 
of non-equilibrium thermodynamics due to a suitable redefinition of an energy momentum 
tensor of the "dark" component that respects a local energy conservation. For a recent 
review on the Black hole entropy on scalar-tensor and f{R) gravity, see Ref. 49|. 

In general relativity, the Bekenstein-Hawking horizon entropy is expressed as 5* = 
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FIG. 5: Cosmological evolution of ^«DE in Eq. (I3.12p as a function of the redshift z for /3 = 1.1 
(dashed hne), /3 = 1.8 (thick sohd hne) and /? = 2.5 (dashed and single-dotted line), where the 
thin solid line shows wde = — 1 (cosmological constant). 



y4/(4G), where A is the area of the apparent horizon 5Ch52|. The Bekenstein-Hawking 

entropy 

A 

(4.1) 



^"4G 



is a global geometric quantity which is proportional to the horizon area A with a constant 
coefficient 1/ (4G). This quantity is not directly affected by the difference of gravitational 
theories. We regard the horizon entropy S in Eq. (14. ip as the one in the equilibrium de- 



scription 



481]. On the other hand, in the context of modified gravity theories including 



/(i?) grav ity a horizon entropy S associated with a Noether charge has been proposed by 
Wald [53|. The Wald entropy 5* is a local quantity defined in terms of quantities on the 
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FIG. 6: Cosniological evolution of Wcs in Eq. ()3.14l) as a function of the redshift z for /? = 1.1 
(dashed hne), /3 = 1.8 (thick sohd line) and (3 = 2.5 (dashed and single-dotted line). 



bifurcate Killing horizon. More specifically, it depends on the variation of the Lagrangian 
density of gravitational theories with respect to the Riemann tensor. This is equivalent to 
S = Aj (4Gcff), where G^s = G/F is the effective gravitational coupling in f{R) gravity 54 1. 
Therefore, we use the Wald entropy in the exponential gravity in Eq. 02.31) 

(1 _ /3e-^/^=) A 



S 



4G 



(4.2) 



In what follows, a hat denotes the quantity in the non-equilibrium description of thermody- 
namics. 

It can be shown that the horizon entropy S in the equilibrium description has the following 



13 



relation with S in the non-equihbrium description 48| : 



with 



dS = — ^-p^dS + — 57^ -diS, (4.3) 



,,,._|l^£.-«/«=f. (4,) 



where diS is the new term which can be interpreted as a term of entropy produced in the 
non-equihbrium thermodynamics. The difference between S and S appears in f{R) gravity 
due to dF ^ 0. Note that S is identical to S in general relativity due to F = 1. From 
Eq. (14. 3p . we see that the change of the horizon entropy 5* in the equilibrium framework 
involves the information of both dS and diS in the non-equilibrium framework. 

In Figs. 7, 8 and 9, we show the cosmological evolution of the horizon entropy 5" in 
Eq. (14. 2 p in the non-equilibrium description of thermodynamics and S in Eq. (14. ip in the 
equilibrium description of thermodynamics for /3 = 1.1, 1.8 and 2.5, respectively. In these 
figures, we illustrate the normalized quantities S = S/Sq and S = S/Sq with 5*0 = n/ {GH^) 
being the present value of the horizon entropy S. Furthermore, we also depict the evolution 
oi H = H/Hq. We note that as S* oc H~^, S increases with time as long as H continues to 
decreases to the de Sitter point, in which H becomes a constant. 

In the high z regime {z > 1), since the deviation of the exponential gravity from the 
ACDM model, i.e., general relativity, is very small, the evolution of 5* is similar to that of S. 
In other words, for the high z regime (the higher curvature regime) F(R) = 1 — (3e~^^^^ ~ 1 
because R/Rs ^1. As z decreases (and R also decreases), the deviation of the exponential 
gravity from the ACDM model emerges, i.e., F{R) < 1 and F{R) decreases. Hence, there 
appears a difference between the evolution of S and that of S. Note that S* > S* oc F{R). 
The present values of S are S{z = 0) = 0.90, 0.96 and 0.99 for /3 = 1.1, 1.8 and 2.5, 
respectively. It is clear from Figs. 7, 8 and 9 that both S and S globally increases with time 
for any values of /3. This confirms that the second law of thermodynamics on the apparent 



horizon always holds. 



Starobinsky's model 48 
of S is to that of S. 



'he similar behaviors for both S and S have been obtained in the 
. Furthermore, we see that the larger /3 is, the closer the evolution 
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FIG. 7: Cosmological evolutions of 5 = SjS^ (solid line), S = S/S'o (dashed line) and H = H/Hq 
(dashed and single-dotted line) as functions of the redshift z for /3 = 1.1. 

V. CONCLUSIONS 

In the present paper, we have studied the cosmological evolution in the exponential 
gravity. We have summarized various viability conditions and explicitly illustrated that 
the late-time cosmic acceleration can be realized after the matter-dominated stage. We 
have also shown that the crossing of the phantom divide from the phantom phase to the 
non-phantom one can occur and the cosmological horizon entropy globally increases with 
time. Phenomenologically, at least in the light of the background cosmological evolution, 
the exponential gravity can be regarded as one of the most promising viable modified grav- 
itational theories because (a) it satisfies all conditions for the viability; (b) in substance it 
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FIG. 8: Legend is the same as Fig. 7 but for (3 = 1.8. 

has only one model parameter; and (c) both the current cosmic acceleration following the 
matter-dominated stage and the crossing of the phantom divide can be realized. 

AcknoAvledgment s 

We thank Mr. Hayato Motohashi for his important suggestions. We are grateful to 
Professor Shinji Tsujikawa and Dr. Antonio De Felice for helpful discussions. The work is 
supported in part by the National Science Council of R.O.C. under: Grant #s: NSC-95- 
2112-M-007-059-MY3 and NSC-98-2112-M-007-008-MY3 and National Tsing Hua University 



16 



t.o 


P= 


1 

=2.5 


1 




1 






1 


1 1 1 


4 


















t 

i 


3.5 


- 
















♦ 

/ 
/ 


3 


















/h 


2.5 


















* 
f 

/ 


2 




















1.5 


. s 


■"^v 






.X 




.' 


.'' 


- 


1 


s 


>> 


^ 




^ 








- 


0.5 


- 


1 


1 




\ 

1 


V 


>s 


1 


1 ' ' "^^"^^ 



-1 -0.5 



0.5 



1.5 



2.5 
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